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Given the general gaussian approximation of (9), the correlation radii take the following 
model independent form for azimuthally symmetric systems [4,6] 



LCMS radii can be found by setting (3l = in (11) and using in eq. (12) a form of 
S(x } K) which is valid in the LCMS frame. We can see that in this frame, R\ L provides 
information about any z-t and/or z-x correlations of the source (as seen in the LCMS). 
Obviously, for the model shown in figure 4, these correlations are not insignificant. 

It has been shown that when the mixing parameters are small, the correlation radii 
simply measure the lengths of homogeneity of the source [9]. For example, for the 
source defined by eqs. (6) and (7), the "side" radius takes the form [4,6] 



Here we can see explicitly that when transverse flow is present (y ^ 0), the length of 
homogeneity measured by the "side" radius is smaller than the geometrical transverse 
radius. Due to the m t dependence, this reduction effect is more pronounced for higher 
momentum and/or larger mass particles. 

Using two specific model sources and a model independent formalism, we have shown 
that there is no a priori reason why an R\ L cross term should be excluded from gaus- 
sian fits to experimental correlation data. Not only will the new parameter reveal more 
information about the source, its inclusion will undoubtedly increase the accuracy of the 
other fitted radii. 
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where we use the notation 






(13) 



3 



where p = \J x 2 + y 2 and v <C 1 is the transverse flow velocity of the fluid at p = Rcj. 

To generate an "out-longitudinal" correlation function, we plug in some numbers and 
then numerically integrate eq. (3), using (6) and (7). For simplicity, we consider a pion 
source with no transverse flow (y = 0) which freezes out instantaneously (St = 0) with the 
following other source parameters: Rq = 3 fm, r = 4 fm/c, 8r] = 1.5, and T = 150 MeV. 
In [6] we show that when the emission function (6) with these parameters is integrated 
over spacetime, it produces a very reasonable one-particle distribution. Figure 3 shows the 
"out-longitudinal" correlator for pairs with Y = —2, K± = 200 MeV, and q s = 0. Figure 
4 shows the same correlation calculated in the LCMS frame, which is the longitudinally 
boosted frame defined by (3l = [2]. Since both figures feature rotations of the major 
and minor axes, it is apparent that the correlations both in the fixed and LCMS frames 
would be much better fit by eq. (2) than by eq. (1). 



400 



120 



200 
q L 
-200 
-400 

-120 -60 60 120 

Pout ( MeV ) 

Figure 3: Correlation calculated numeri- 
cally using (3), (6) and (7) in the fixed 
center of mass frame. 
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Figure 4: Same as fig. 3, but calculated in 
the LCMS frame. 



4. Model Independent Derivation of the Correlation Function 

By making a saddle point approximation to a general emission function, we can de- 
rive model independent expressions for all of the correlation radii [4,6]. We define the 
spacetime saddle point x of the emission function S(x } K) through the four equations 

\nS(x,K) 



dx 
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where p = {0,1,2,3}. Essentially the saddle point is that point in spacetime which 
has the maximum probability of emitting a particle with momentum K. A saddle point 
approximation for S(x } K) can then be made in the following way [6] 



S(x, K) ~ S(x, K) exp 



E 



[X, 



x n) 



2A^ 



where we define lengths of homogeneity and source mixing parameters by 



K(K) 



In S(x, K) 



-1/2 



and B^(K) 



d d 



In S(x, K) 



(9) 



(10) 
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2. Stationary Gaussain Sources 



Consider the following azimuthally symmetric gaussian emission function: 

z 2 (t-t ) 2 ' 



S(x, K) = f(K) exp 



2 i 2 

x + y 



(4) 



2R 2 2L 2 2(St) 2 

Using (3), one can see that the corresponding correlation function takes the form 

C(q, K) = 1 ± exp[-q 2 R 2 - q 2 L (R 2 + (3l(St) 2 ) - q 2 L (L 2 + (3 2 L (Stf) - 2q L q L (3 L (3 L (St) 2 ] , (5) 

so the q±qL cross term provides a measurement of the duration of particle emission (St). 

The effect of the cross term is easiest to see by looking at the "out-longitudinal" pro- 
jection of the correlation function. Figures 1 and 2 show contour plots of the "out- 
longitudinal" correlator for a source featuring R = L = cSt = 3 fm and pairs characterized 
by q s = 0, j3]_ = \CE c. The outer contours are for \C — 1| = 0.1, and each successively 
smaller contour represents an increase of 0.1. In figure 1/3^ = 0, so the cross term van- 
ishes, thus causing the major and minor axes of the elliptical contours to coincide with 
the "out" and "longitudinal" axes. In figure 2, however, (3^ = \CE c, and the nonvanishing 
cross term causes a rotation in the major and minor axes. This effect should be easily 
observable, and has in fact already been measured in NA35 correlation data [7]. 
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Figure 1: Eq. (5) with (3 L = 0. Figure 2: Eq. (5) with (3 L = ^j\j2c. 

3. Model Featuring Boost-Invariant Expansion 

Now we turn to a more realistic model, similar to those in [8]. In the center of mass 
frame of an expanding fireball, we define the following emission function [4,6]: 



S(x, K 



T m t c\i(ri — Y) 



(2it) 3 tJ2it(St 



exp 



K-u(x) 



x 2 + y 2 
~2RY 



TO 



2(Sr]) 2 2(StY 



(6) 



where T is a constant freeze-out temperature, r = \/t 2 — z 2 , rj = | \n[(t -\- z) / (t — z)], m t = 

\J m 2 + Kj_ , and Y is the rapidity of a particle with momentum K. We consider a flow 
which is non-relativistic transversally but which exhibits Bjorken expansion longitudinally, 

u(x) ~ ( (l + \(vplR G f) eh*/, (vx/Rg), (vy/Ra), (l + \(vplR G f) sbrA , (7) 
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q Using two specific models and a model independent formalism, we show that an "out- 
longitudinal" cross term should be included in any gaussian fits to correlation data. In 
addition, we show that correlation radii (including the cross term) measure lengths of 

^ homogeneity within the source, not necessarily geometric sizes. 

q 1. Introduction 
IT) 

Experimentally measured Hanbury-Brown Twiss (HBT) correlations between two iden- 
tical particles are typically fit by gaussians of the form [f ,2] 

^C( q ,K) = l±Xe X p[-q 2 s R 2 s (K)-qlRl(K)-q 2 L R 2 L (K)} , (f) 

where q = pi — p 2 , K = |(pi + P2), the + ( — ) sign is for bosons (fermions), and the 
coordinate system is defined as follows [3]: The "longitudinal" or z (subscript L) direction 
is parallel to the beam; the "out" or x (subscript _L) direction is parallel to the component 
of K which is perpendicular to the beam; and the "side" or y (subscript s) direction is 
the remaining direction. We claim that significantly more can be learned and better fits 
achieved if an "out-longitudinal" cross term is included in the following way [4] 

C(q,K) = i ± Aexp [-q 2 s R 2 s (K) - qlR\{K) - q 2 L R 2 L (K) - 2q L q L Rl L (K] 

where R\ L is a parameter which can be either positive or negative. 

To see how this cross term arises in two-particle correlations, we use the following well 
established theoretical approximation [5,6] 

C(q ,K)c 1± '^;^;;;' 2 , < 3 > 

I J d 4 x b(x, R)\ 2 



where q = E\ — E 2 and K = Ek = \Jm 2 + |K| 2 . Here S(x,K) is a function which 
describes the phase space density of the emitting source. For pairs with |q| <^ Ek, we 
can use the approximation q — (3±q± -\- ^lIl, where = Ki/ Ek- 



